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I. INTRODUCTION 

In the Standard Model, CP violation is predicted in 
weak decays thanks to the single phase of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix. It is also well known 
that such a weak phase is not sufficient to generate a CP 
violating decay amplitude. Strong phases are necessary 
and their strength may significantly enhance the effect 
of the weak phase. Therefore, hadronic effects, such as 
resonances of daughter particles in 5*- and higher waves, 
require a careful analysis in the determination of CP vio- 
lating phases in hadronic two- and three-body decays [1- 
4]- 

The antimatter-matter asymmetry is expected to be 
very small in weak decays of B s mesons; any observed 
deviation may well be a signal of physics whose origins 
lie beyond the Standard Model. In the B® — > J/ip<j> chan- 
nel, recent measurements by the CDF [5] and D0 [6, 7] 
Collaborations of the B s — B s mixing phase, — 2/3 s , while 
not definitive, are considerably larger than Standard 
Model predictions. Taking advantage of the fact that 
the B° s -> J/V>/ (980) channel does not require any an- 
gular analysis, one can compute the ratio between the 
B° -> J/ip(j) and B° s -> J/V>/ (980) decay widths in or- 
der to estimate the tt + tt~ S'-wave effect on the value of 
j3 s . A first qualitative attempt to predict the ratio, 

T(B ( > -» J/^/q(980),/o(980) -» tt+tt-) 

/o/0 r(B»4j/#,^M-) ' ( ' 

was made by Stone and Zhang [8] and gives a result of 
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the order of 20% — 30%. Their estimate relies on exper- 
imental data on D+ -> / (980)7r+ and D+ -> </>7r+ de- 
cays and seems to indicate that the 5-wave contribution 
of /o(980) K + K~ cannot be ignored when analyzing 
the angle (3 S in B® — > J/ip<j>. Likewise, Xie et al. found 
the effect of an 5*-wave component on 2(3 S to be of the 
order of 10% in the <j> resonance region [9]. 

Based on the QCD factorization (QCDF) formalism we 
perform a first robust calculation of the ratio 1Zf Q i^. To 
this end, all the available obscrvablcs (polarizations and 
branching ratio in B® — > J/ip<t>) are used to effectively 
constrain the analysis of the B® — > J /ip<t> channel. The 
branching ratio and CP asymmetry are then predicted 
for Bg — > J/-0/o(98O), where we assume that merely the 
ss component of the /o(980) is involved in the hadronic 
B s — > / (980) transition matrix element. 

In Section II we introduce the general expressions for 
the B° -> J/ipcj) and B° -> J/-0/ o (98O) weak decay am- 
plitudes whereas Sections III and IV provide the details 
on the leading order corrections in a s for both these 
amplitudes, respectively. In Section V, we list all nu- 
merical values of input parameters and briefly recall our 
model for the B s -4- /o(980) transition form factor [10] 
on which the ratio Ttf /tf> directly depends; we also define 
the paramctrization for the B s — > cp form factor. Sec- 
tion VI is devoted to our results and, finally, conclusions 
are drawn in Section VII. 



II. GENERAL FORM OF THE B° cbj/ip AND 
B° ->■ /o(980)J/V DECAY AMPLITUDES 

It is important to realize beforehand that the applica- 
tion of QCDF, following Refs. [11-14], to B° s decays into 
a heavy-light final state is not self-evident. In both fi- 
nal states, (j)J/tj) and /o(980) J/tp, the s-spectator quark 
is absorbed by the light meson while the emitted me- 
son is heavy in which case QCDF is not reliable [11]. 
Nonetheless, as argued in Refs. [15, 16] and more re- 



2 



cently in Ref. [17] , the production of a heavy charmonium 
qq pair bears "color transparency" properties similar to 
those of a light meson, provided this color-singlet pair is 
small compared to the inverse strong interaction scale, 
1/Aqcd- This was explicitly demonstrated in next-to- 
lcading order calculations for exclusive B decays to J/tp 
final states (J/ipK, J/ipK*), where infrared divergences 



were shown to cancel [15, 16]. 

In the following, we present the B® decay amplitudes in 
which the short- and long-distance contributions are fac- 
torized in the approximation of a quasi two-body state, 
MiM 2 , where either M 1 M 2 = fo (980) J/ V> or M 1 M 2 = 
(f>J/ip. We begin with the B° s — >■ <pJ/tp amplitude which 
can be written for each helicity, h = —1,0, 1, as [14], 



Ah 



h 

4>J/i> 



q—u,c 



(a q 2 h (m b ) + C 71 ) + 4 h (m b ) + af h (m b ) + af l (m b ) + a q 9 h (m b ) 



Summing over all the possible helicities, the squared modulus of the total amplitude reads 



A 



B°-Ht>J/il>\ 



\A 



h=-l 



i I Ah=0 
"I- \-^B°^J/^\ 



+ \A h=+1 I 

\-^B^<pJ/tp\ 



(2) 



(3) 



The B® — > <j)J/ij} decay amplitude is obtained by exchange of helicity signs, h = +1 — > h = —1, and replacing X q by 
its complex conjugate. The — > /o(980)J/^> amplitude is, 



■A B o^f j/4, = X i\ A/W/V> V( a 2( m b)+C) +a q 3 (m b ) + a|(m fc ) + a q 7 {m b ) + ag(m 6 ) 



(4) 



foJ/i> 



The different elements entering in the amplitudes (2) and 
(4) are defined in Eqs. (6), (7), (15), (22) and (24). The 
CP conjugate B° s decay amplitude is again found by re- 
placing X q by its complex conjugate. 

With the generic amplitude, Ab^^m^j/^, the branch- 
ing ratio, 



B(B° S -> M x Jji>) = — — 

v s ' ' T B o 167rm B o 

x A 1 / 2 ^1 , m 2 Mi I m 2 B o , m 2 j/^ / m 2 B o ) \A B o^ MlJ/i> \ 2 , (5) 



can be computed. The J/ip mass is noted mj/^ while 
mM, = w/ (980) or m 4> denote the /o(980) and 4> masses; 
the triangle function is X(x, y,z) = (x + y — z) 2 — 4xy. 
In Eq. (5), r B o = 1/t b o is the B® decay width with 
t b o = (1.470 ± 8 0.026) ps S [18] and m B o is the B° s mass. 
For the CKM elements in Eqs. (2) and (4) we use the 
Wolfenstein parametrization, 



X u 

X r . 



V: b V us = AX A {p + ir 1 ) , 
A 2 



v: b v cs 



AX 



(6) 



with the Wolfenstein parameters A = 0.814, p = 0.1385, 
r] = 0.358 and A = 0.2257 [18]. 



A. Non-perturbative amplitude 

1. The case of the scalar-vector decay 
The scalar-vector factor, Af j/$, in Eq. (4) is given by, 

A foJ/4 , = (/ob/ )lH(l ~ J 5 )s\B° s ( PB o)) 

x <J/^(p w ,e^)|cyc|0) , (7) 

where the hadronic matrix element which describes the 
transition between the B® and a scalar meson, / , with 
the respective four- momenta p B o and p/ is [19], 



(f (p fo )\b ltl (l- 75 )s\B s (p B o)) 



(pso 



m B o — m 



^0) if °^V) 



m 2 R0 - m 2 f ,,,, ,. 
+ \ 2 N m^V), (8) 



with q = p B o - pf Q , q 2 = m 2 }/ ^ and where F ± s °(q 2 ) 

B°— y f 

and Fg ' J0 {q 2 ) are the vector and scalar form factors, 
respectively. In Eq. (7), the leptonic decay constant, 
/ju, of the J/ip vector meson, with four-momentum, 
Pj/ip, and polarisation, e*j^, is defined as, 

(•//V'(Pj/^,£*//^)|c7 ai c|0> = -ifj/^mj/^*^ . (9) 

The scalar-vector factor, given by the product of 
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Eqs. (8) and (9), is then obtained as, 

A f J/->p = - i —^ 2m J/4> e j/</> ' Pb« F 1 s {mjfy)fj/4, , 

(10) 

with 4m^|e^ • p B °\ 2 = m 2 B o X 1/2 (m 2 B o,m 2 j^,m 2 fo ) 
and the Fermi constant, G F = 1.16 x 1CT 5 GcV~ 2 . The 
B° — > /o transition form factor s °( TO j/^) wm be 
discussed in Section V. 

2. The case of the vector-vector decay 

For the case of two vector mesons, Mi and M 2 , the 
helicity formalism requires the introduction of three po- 
larization four-vectors, eu^.k (j = 1,2 and k = 1, 2, 3) for 
each spin-1 particle, Mj, 

e Mj,i = (0, e*Mj,i) , 

e Mj,2 — (0, ^Mj,l) i 

£M„3 = {{PMjl/mM^EMjPMj/mMj) ■ (11) 

where tom.,- , PMj and i?Mj are the mass, the momentum 
and the energy of the vector meson, Mj, respectively. 
The energies E Ml , E M2 are given by, 

e m 1<2 = { m B° ~ m lh ~ m k) • ( 12 ) 

In Eq. (11), pMj is defined as the unit vector along the 
momentum: p Mj =PM j /\pM j \- 



The three polarization four- vectors, eMj,k, also satisfy 
the following relations, 

eM 3 ,fc 2 = -1 , and e Mj ,k ■ e Mj ,i = , for k ^ I . (13) 

The vectors cm^i, ^Mj,2 and cm 3 -, 3 form an orthogonal 
basis in which £Mj,i and 6^,2 describe the transverse 
polarizations while e*M 3 J 3 is the longitudinal polarization 
vector. With these three vectors one builds up the helic- 
ity basis, 

= ^ (eMj.i +* %,2) = ^(0,+l,i,0) , 
£ m 3 - = ^= ( £ M j: i - i ^,2) = -^=(0,-1, i,0) , 

£ Mj,0 = EMj,3 ■ (14) 

and £mi,± = £ m 2 ,=f- ^ n Eq. (14), the new four- vectors 
and eMj,o are eigenvectors of the helicity 
operator corresponding to the eigenvalues h = +1,-1 
and 0, respectively. 

The vector-vector factor, A MiM2 , in Eq. (2) is 

A H Mi m 2 = {Mi{PM 1 ,e* Ml )\b ltl {l-^)q\B Q s {p B o)) 

x (M 2 (p M2 , £^ 2 )|f7V|0) , (15) 
where, in the B® rest-frame, the vector mesons Mi and 
M 2 have opposite momentum p Ml = —j>m 2 along the 
^-direction and cm.,-,0 ■ PMj = 0. 

The matrix hadronic element of a P — > V transition 
can be decomposed into Lorentz invariants as [16, 19, 20] 



A B °~* Mi ( 2 ) 

(M j (p M j,e* Mj )\b'y ll (l-'Y5)q\B°(pBo)) = e* M .^(m B o + m Mj )A Bs ^ Mj (q 2 ) - (p B o +Pm j )^(s* Mj ■ p B o) ^— -j-^— 



1 ■ *v ql p 



2V 



(q 2 ) 



m B o + m Mi 



(16) 



r 



where the form factors Aq s ~^ M ' (q 2 ), A Bs ^ Mj (q 2 ), 



,B"^Mi 



A^"^ 1 ™ 3 (q 2 ) and ^4 3 s ^ 3 (q 2 ) obey the following exact 
relations, 



A?^{q 2 ) 



™^z™^ A ?^ [e) , (17) 

Zm Mi 



as well as for q 2 = 0, A B °°^ Mi (0) = A B °^ M *{0). 

Specifically for Mi = <fr, and M 2 = J/ip, the helicity 
dependent vector-vector factor A^j^ in Eq. (2) has thus 



the following form, 



A (h=0) _ i^F fl/l 



-m^mBO+m^A^ {m 2 J/4 ,) 



+ (m|o + ml - m^/^Af^^im 2 ,^) 



(18a) 



Ah=±i) 



"*JH> - ^ m B°mj/^f.J/4,F T s (m J/v ,) . (18b) 
In Eq. (18b), the transition form factors F± B ^^(q 2 = 



4 



(™W = I 1 + ) ( m JA/>) 



m,B° 



m B + 777-0 



^^K/V-) , (19) 



where the center-of-mass momentum \pb°\ is defined as, 



Wb° I - 



(m|o - M|) (m|o - Mi) 

2777 B 



(20) 



with M± = mju ± 777,0. We note that a somewhat dif- 
ferent form for A^Jj^j was derived in Ref. [20], which 
seems to approximate the vector mesons as light mesons. 
The form factors Aq 3 ~ > (trij^) and A 1 "^ f (rrij,^) in 



Eqs. (18a) and (19), as well as V B * 

are defined in Section V. Ref. [14] asserts that when 
neglecting vector meson masses, Eq. (18a) reduces to, 



* ^/j/^™|oA^ (777 2 7 ) . (21) 



y/2- 



The numerical effects in the calculated values of B® —> 
J/ip(f> and Bg — > J/t/>/ (980) branching ratios are too 
important to justify such an approximation. 



B. Perturbative amplitude 

The a^ h (fi) coefficients that appear in Eqs. (2) and (4) 
are linear combinations of Wilson coefficients, C„(/x), ei- 
ther at the scale /x = mj, or mi,/ 2 (see below): 



a* h (m b ) 



C n (m b ) + 



C n ±i(m b ) 



N r . 



N n (J/iP) 



47T 



+ rrC F a s ( mb /2) Cn±l{ ™ b/2) H^J/iP) . (22) 

The superscript, (h), explicits the helicity dependence 
of a^ h (iJ,) in the case where B® decays into two vector 
mesons. This superscript is dropped in the scalar-vector 
case. There is no flavor dependence in a^ h (fj,) for n = 
1,2. In Eq. (22), the upper (lower) signs in C n ±i(fi) 
apply when 77 is odd (even) and 

N n {J/ij) = 0, 77 e {6, 8}, else N n (J/ip) = 1 . (23) 

The Wilson coefficients, C„(/i), in the Naive Dimensional 
Regularization (NDR) scheme are taken at the hard scale 
777;, for the vertex, V^(J/ip), and penguin, P%' h (J/i/j), 
corrections, whereas in the hard scattering, H^MxJ/ij)), 
amplitudes they are evaluated at m b /2 since those con- 
tributions involve the spectator quark. The strong cou- 
pling constants at these scales are a s (m b ) = 0.224 and 



a s (m b /2) = 0.286 [18], while the number of active fla- 
vors is rip = 5, the color number iV c = 3 and Cp — 
[Nl - 1)/2N C . 



C. Suppressed higher order corrections and 
possibility of new physics 

There are no contributions, such as given by the anni- 
hilation operators derived in Ref. [13], to the two decays 
considered here. This is because for the final states, J/ip(j> 
and /o (980) , both mesons are simultaneously flavor 
and color singlets. At tree level, for instance, the W ex- 
change diagram produces the charmonium cc, yet the cre- 
ation of the ss which hadronizes to an /o(980) or <fi must 
proceed via multiple gluons or by means of photon/Z ex- 
change. The annihilation is thus either strongly (Zweig) 
suppressed in a s or the suppression is in the electromag- 
netic coupling constant a om . 

On the other hand, as will be discussed in Section VI, 
if we account for vertex, penguin and hard scattering 
corrections only, the B° s — > J/ip(f> observables are only 
moderately well reproduced. As can be seen in Table IX, 
the branching ratio, for instance, is about 20% too large 
(although still within the experimental errors). We there- 
fore allow for additional phcnomcnological amplitudes 
that mock up "other" contributions, be it from annihi- 
lation topologies expected to be strongly suppressed or 
possible physics beyond the Standard Model [21]. These 
arc included in Eqs. (2) and (4) with the amplitudes, ( h 
and (, conveniently scaled as, 



Bmi J/ip 
^Mi J/ip 



X c ■ 



(24) 



The factor B MlJ /^ is chosen to be a product of decay 
constants, either 



G F 

B f j/1> = ~i /bo ff fj/i, 



if Mi = / (980) or 



B <t,j/^ - i /bo 



,/. 



J/V 



(25) 



(26) 



if Mi — <fi, while the factor Xc is a complex parame- 
ter discussed in Section VC. We note that the decay 
constant, ff , vanishes due to charge conjugation invari- 
ancc, wherefore the scalar light cone distributions ampli- 
tude (LCDA) is normalized to ff — ff rrif /(m Uy d(^) — 
Wtt,d(/j)), which is finite [22]. We shall return to this 
issue in Section IV. 



D. The ratio 7l fo /^, 

Prior to discussing the various a s (/i) corrections to the 
amplitudes, a^ h {fj), it may be of interest to observe the 
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qualitative behavior of the ratio, TZf / ( j > , in terms of the 
scales Aqcd and mj. A naive factorization analysis yields 
a hierarchy of helicity amplitudes for B into vector- vector 
decays [14], 



A (h=o) 



. A {h=+l) 

■ ^BO^tjiJ/ip ■ • n B°-t<l>J/il> 



1 : 



• A (h= ~ 1) 

Aqcd 
m b 



Aqcd V 
m b ) 



(27) 



while for B s mesons the signs are exchanged (h = +1 



(h=0) 



and 



h = —1). Furthermore, the amplitudes A 

"4b°->/ •//'!/> are °f same order in Aqcd/"^- With this 
estimation, the ratio 7£/ o /0 we are interested in becomes, 



two particle distribution amplitudes are approximated by 
4> v {x) =6x(l-x) (29) 

and 

<p v (x) = 3(2x - 1) , (30) 

respectively. In the annihilation and hard scattering am- 
plitudes, the chiral coefficient, r^(/x), is defined as 



Vi \ _ 2m v fy(p) ^ 2m v 
x rn b (fi) f v ~ m b (n) 



(31) 



11 



fo/4> ~ 



>JV/V> 



(h=0) 



B°- 



p , utfc-i) l 2 ^U(' l=+1 ) I 2 

4>J/il>\ " r l^BO-HfrJ/ip] " r \-"-B°^r<t>J/il>\ 



0(l)+0 



A 



QCD 

m b 



o 



A, 



QCD 

m b 



(28) 



Hence, -R/ o /0 is 0(1) for Aqcd/w;, corrections. 

Nonetheless, non-perturbative hadronic effects can 
spoil the naive factorization and violate the hierarchy 
in Eq. (27); so do electromagnetic penguin contributions 
where a photon with small virtuality subsequently con- 
verts into a vector meson [23]. 



III. QCDF CORRECTIONS FOR B° cfrJ/ip 
DECAY AMPLITUDES 

Due to the structure of the four-quark operators in 
heavy quark effective theory and the conservation of the 
flavor quantum numbers, the final state M X M 2 = <pJ/"4> is 
created from the transition B® — > (f> and the production of 
J/ip from vacuum. As discussed in Section II, the decay 
amplitudes at leading order in K-Qcn/mb and a s (m b ) are 
given by the factorized product of a transition form factor 
and a decay constant. Following Rcf. [14], we only give 
QCD corrections that explicitly appear in the amplitude 



Ah 



of Eq. (2). 



We discard terms proportional to r — (mj/^/mBj 2 — 
1/3 in vertex corrections which stem from the presence 
of the charm quark in the loop diagram; we have nu- 
merically checked that their contributions to the a*' ft (/x) 
coefficients are negligible, all the more so when seen in 
the light of the large hadronic uncertainties of the form 
factors [see Sections (VA) and (VB)]. We note that in 
the limit r — > 0, one recovers the vertex correction known 
from, for example, B — > irn which is of course infrared 
safe. 

Since the coefficients in the Gegenbauer expansion of 
the LCDA are poorly known for the scalar mesons, and 
only with non- negligible errors for the vector mesons V — 
(j) and V = J '/ip, we limit ourselves to leading terms in the 
expansion. The leading twist-2 distribution and twist-3 



where fy (/z) is the transverse decay constant for any vec- 
tor V and /z = m b /2. 



A. Penguin contributions 

The penguin contributions to the amplitude in Eq. (2) 

stems from the positive helicity, h = +1, amplitudes 
p 9 ,/i=+i 

^7,9 



(J/VO given in Ref. [14], 



m B om b 2a e 

H 

m%,, 27n 



(C^) + N c C 2 (n)) 



J/V 

2 2 
Tfl V 

S qc In — f + Sgu In + 1 
zr fj, 2 



whereas P^'g ""V/V) = 0. In Eq. (32), fj, 
G* (n) = CV 7 (/x) - C 5 (zt)/3 - C 6 (zt), a e = 1/129 is 
the electromagnetic coupling constant and the scale v 
refers to the /jm decay constant scale. One also has 



(32) 



3,5 



(J/V>) = as well as Pl'j h 7 ~g\j/^) = 



Vertex contributions 



In B® — > (f>J/ip, the electroweak vertex receives a s (fi) 
corrections to all a^ h (^i) in the amplitudes Ago^j^. 
For h = 0, these are, with /z = m b , 



12 In 



mi, 



'din 



27 



KTV/v) = 



zz j 2 ' 

for n e {2, 3, 9} 

for n G {5, 7} 



(33) 
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whereas for h — — 1 one has, 



12 In 



143 



M 7 4 
for n e {2, 3, 9} 



-12 In (?* 



7T 2 + 



95 



for n e {5, 7} 



and for /i = +1 one has, 



(34) 



for n e {2, 3, 9} 



-12 1n^-j- Y+ 6 i7 r +T , 
for n e {5, 7} . 



(35) 



where B 1 (m b /2) = -0.54 [22] is the / (980)'s first Gegen- 
bauer moment and we remind that only odd moments 
contribute in case of charge-neutral scalar mesons. In 
particular, contrary to the pseudoscalar LCDA, the lead- 
ing term 6x(l — x)B vanishes since B = (mi(/i) — 
™>2{lJ)) / m s, where ms is the scalar meson mass and 
™i.2(m) its running quark masses. The scalar twist-3 
two-particle distribution is given by 



<ff {x) = 1 



(40) 



The asymptotic forms of the LCDA, (pj/^{x) (Eq. (29)) 
and <pj/^(x) (Eq. (30)), are used. As in the B® — > (pJ/tp 
decay, the J/tp meson is created from vacuum whereas 
the transition B® — > / (980) produces the scalar meson. 
Here, we only consider the ss component of the /o(980) 
since the flavor of the spectator quark in the tree and 
penguin topologies of B® decays is strange. There are no 
penguin corrections [13] to the B® — y /o (980) J/V> decay 
amplitude in Eq. (4). 



C. Hard scattering contributions 



Vertex contributions 



The gluon exchange between a J/tp meson and the 
spectator s-quark leads to the hard scattering ampli- 
tudes, 

H'r°(^JM = ± 3^^(r*MX B + 3), (36) 

for h = 0, fi = TO h /2 and A B o = 0.350 GcV [13]. The plus 
sign is for n = 2, 3, 9 and the minus sign for n = 5,7. The 
phcnomenological amplitude, Xh, parametrizes the end- 
point divergence of the scalar meson's LCDA and is de- 
fined in Eq. (54). For the helicity, h = +1, the correction 
reads 

i/r +1 (W)-Tl8§^^^(A ff -l), (37) 

where the minus sign applies to n — 2,3,9 and the plus 
sign to n = 5,7. The helicity, h = —1, contribution is 
simply, 

H%=- 1 (<j)J/il)) = Q for n = 2,3,5,7,9 . (38) 



IV. QCDF CORRECTIONS FOR 
B° ->■ /o(980)J/V DECAY AMPLITUDES 

We now turn to the B® — > J/tpfo(980) transition 
for which the a s (fi) corrections are all included follow- 
ing Ref. [13] applied to an SV final state. For previ- 
ously mentioned reasons, we solely employ the first non- 
vanishing leading term in the LCDA, 



(j) So {x) = 6x(l - x) 3B 1 (p)(2x - 1) , 



(39) 



At the order of a s (/z), the vertex correction, V n (J/tp), 
involves the leading twist distribution, (pj/^,(x), and a 
gluon kernel given in [13]. We derive from this the ex- 
pressions, 



v n (JH>) = { 



for n e {2, 3, 9} 



(41) 



for n E {5, 7} 



with ji = m b . 



B. Hard scattering contributions 

The hard scattering correction in case of an foJ/tp final 
state reads 



ff„(/ J/V) = ±3 



B f J/i> m B» f- fo 

A f J/i> 



(42) 

where the plus sign applies to n = 2, 3, 9, the minus sign 
to n = 5, 7 and Xh is given, as in the case of the (pJ/tp 
final state, by Eq. (54). 

The chiral coefficient, r£°(/x), enters Eq. (42) rather 
than r£° (//) defined as, 



2m 2 , 

JO 



m b (n) (mi(/z) - m 2 (/i)) 



(43) 
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TABLE I: Wilson coefficients at the fj, = mb and p — mt/2 scales in the NDR scheme [12]. The coefficients CV(/i) — Cio(p) 
must be multiplied by a e . 





din) C 2 (p) Cs(Ai) C 4 (A«) C 5 (A*) C 6 (/i) CV(a«) C 8 (p) C 9 (/i) Cio(aO CV 7 (/i) 


p = m b 
p = m b /2 


1.081 -0.190 0.014 -0.036 0.009 -0.042 -0.011 0.06 -1.254 0.233 -0.318 
1.137 -0.295 0.021 -0.051 0.010 -0.065 -0.24 0.096 -1.325 0.331 -0.364 



TABLE II: Values of the higher order correction (pc, <t>c) 
and hard-scattering (pH, 4>h) parameters as function of the 
B° decay constant. 



f B o [MeV] 


Pc 4>c (°) Ph 4>h (°) 


230 
260 
290 


4.52 ±2.24 173.8 ±37.6 1.90 ±0.20 266.0 ±21.6 
6.16 ±2.03 176.1 ±53.6 1.70 ±0.16 260.6 ± 19.3 
7.33 ± 1.63 176.0 ±57.6 1.54 ±0.15 255.6 ±17.3 



The reason is that in case of neutral scalar mesons, 
mi(/x) = m2(/i) and r£°(/i) diverges. On the other hand, 
it is known from C-conjugation invariance that the vector 
decay constant of the neutral scalar meson must vanish. 
However, the quark equations of motions yield a relation 
between the scalar- and vector-decay constants, ff and 
ffo respectively: 



ffo ~ 



m fo 



mi(,u) - m 2 (/u) 



ffo 



(44) 



where to/ // = (0|<f2<?i|/o)- Since //„ is non-zero, the 
product f f m / / '(mi(/j) — m2(/x)) is finite in the limit 



with 



m 2 (/i). We thus recombinc, ff r x 



fo 



ffo r 



fo 



rfo 



(p) 



2m fo 

m b {n) 



(45) 



V. NUMERICAL PARAMETERS 

This section serves to summarize all parameter values 
required for numerical applications. The Wilson coeffi- 
cients at the scales p, = m b and fi = m b /2 used in this 
work are listed in Table I. For the meson masses, we refer 
to the latest PDG values [18], which are (in GeV): 

m B o = 5.366 , m B * = 5.412 , m fo = 0.980 , 

mj/1, = 3.096 , to = 1.019 . (46) 

The running quark masses at \i = m b = 4.2 GeV are (in 
GeV), 



and those at n = m b /2 = 2.1 GeV are, 

m b = 4.95 , m c = 1.51 , m s = 0.09 , m u . d = 0.005 . (48) 

We take the <f> decay constant values from Rcf. [14]: 
/ = (221 ± 3) MeV and f£ = (175 ± 25) MeV. For the 
J/iJj meson, we use f J/4 , = (416±6) MeV [24] and fj,^ = 
(405 ± 5) MeV [16]. In the B s J/^/ o (980) channel, 
the ss component of the /o(980) is involved which implies 
the poorly known scalar decay constant //„ : one theoret- 
ical estimate yields fp = (180 ± 15) MeV [25] whereas 
a much larger value "f fs (1 GeV) = (370 ± 20) MeV 
[//•(2.1 GeV) = (460 ± 25)MeV] is found in Rcf. [22], 
both from coupling to the scalar ss current only (de- 
noted by the superscript s in /q, which we use hence- 
forth). Similarly several theoretical predictions exist 
for the leptonic B s decay constants of which we se- 
lect three values from unquenched lattice QCD: f B o = 
(204 ± 12±||) MeV [26], f B o = (259 ± 32) MeV [27] and 
f B o = (231 ± 15) MeV [28]. 3 

To illustrate the sensitivity of the ratio 7£/ m to the 
hadronic uncertainties, we exemplarily choose three dif- 
ferent values for each decay constant: f B o = 230, 260, 
290 MeV and j f . = 340, 380, 420 MeV. 

A. B — ¥ V transition form factor 

Values for the J5° — > 4> transition form factors are taken 
from the pole-extrapolation model by Melikhov [19]: 



A (q 2 ) 



2\B" 



oo(0) 



1 



l-a 1: 



b5 



(49) 

The form factor V(q 2 ) Bs ^^ is given by a similar ex- 
pression in which a (0) is replaced by v(0) and m B o 

by m B * [19]. Next, the A 1 (q 2 ) B °^' f > form factor is 
parametrized by 



2\B" 



oi(0) 



1-<7 1; 



(50) 



Finally, A 2 (q 2 ) B ^ 4 ' has the same functional form as 
Ai(q 2 ) B ^^ where ai(0) is replaced by 02(0). In both, 



m b = 4.2 , m c = 1.3 , m s = 0.07 , m u . d = 0.003 , (47) Eqs. (49) and (50), the momentum transfer is q 2 



rn 
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TABLE III: Short-distance amplitudes, a^' h (m b ) x 10 3 , for the helicity h = +1 in — > J/V>0, as a function of the decay 
constant, f B o, and with //« = 380 MeV. The LOVP results are obtained with the leading order (LO) amplitude to which 
vertex V and penguin P corrections are added. In case of LOVPH, the hard scattering contribution with the endpoint 
parametrization X H is included. LOVPH + C contains additionally the purely phenomenological contribution f C 1 ) with two 
more parameters. 



/so [MeV] 




230 


260 


290 




LOVP 


LOVPH 


LOVPH + C 


LOVPH 


LOVPH + C 


LOVPH 


LOVPH + C 


u.c / \ 

a 2 ( m b) 


60.38 -i 161.7 


-3.77 + i 148.8 


-8.43 + 4 129.87 


38.06 + i 149.7 


8.21 + i 130.04 


75.52 + i 140.4 


22.57 + 4 128.49 


U.C { \ 

a 3 ' (nib) 


5.66 + i 5.39 


8.54 -i 8.54 


8.75 -i 7.69 


6.66-i8.58 


8.0-4'7.70 


4.98-4'8.17 


7.36 - 4 7.63 


U.C / \ 

a 5 ( m b) 


-5.27 -i 6.28 


-8.94 + i 11.47 


-9.21 + i 10.39 


-6.55 + 4-11.52 


-8.25 + 4 10.40 


-4.41 +4 10.99 


-7.43 + 4 10.31 


a%(m b ) 


0.12 + i 0.07 


0.17- *0.13 


0.17-i0.12 


0.14-40.13 


0.16 -i 0.12 


0.11-40.13 


0.15-40.12 


ar{m b ) 


0.69 + i 0.07 


0.73 -i 0.13 


0.74-i0.12 


0.71 -i 0.13 


0.73-i0.12 


0.68 - 4 0.13 


0.72-4'0.12 


ag(rn b ) 


-9.25 -i 0.27 


-9.40 + 4 0.43 


-9.41 +4 0.39 


-9.30 + i0.43 


-9.37 + 4 0.39 


-9.22 + 4 0.41 


-9.34 + 4-0.38 


ag(m b ) 


-8.68 -i 0.27 


-8.83 + 4 0.43 


-8.84 + i0.39 


-8.73 + i 0.43 


-8.80 + 4 0.39 


-8.65 + 4 0.41 


-8.77 + 4-0.38 



TABLE IV: Short-distance amplitudes, a„ h (nib) x 10 3 , for helicity h — —1 and — > J/ip4>. Since the hard-scattering 
contributions are zero, these amplitudes are independent of f B o . 



a2 ,c (m b ) 


-51.72 


a%' c (mb) 


9.39 


U.C / \ 

a 5 ( m b) 


-9.63 


U.C / \ 

a 7 ' (nib) 


0.12 


ag' c (m b ) 


-9.49 



In Eqs. (49) and (50), the form factors at q 2 = 
are a (0) = 0.42 (v(0) = 0.44) and ai(0) = 0.34 
(a 2 (0) = 0.31). The extrapolation parameters are, for 
A (<7 2 ) s °^, cti = 0.55 and o 2 = 0.12; for Viq 2 ) 3 "^, 
a 1 = 0.62 and a 2 = 0.20; for A^q 2 ) 3 "^, cr 1 = 0.73 and 
a 2 = 0.42 and finally for A 2 (q 2 ) B "^' t ', a 1 = 1.30 and <r 2 = 
0.52. The respective values for the form factors at the 
value q 2 = m 2 ^ arc A (q 2 ) B s^ = 0.76, A^q 2 ) 3 "^ = 

0.42, A^q 2 ) 3 "^ = 0.49 and Viq 2 ) 3 "^ = 0.80. 



B. B —¥■ S transition form factor 

We studied the transition form factor, F { (q), 
in a comparative calculation using a dispersion relation 
and a covariant light front dynamics model [10]. To our 
knowledge, this form factor has only been calculated re- 
cently in QCD sum rules [29, 30] and pQCD [31] for q 2 = 



and must be extrapolated to the value F Q (^ (m^). 

In our work [10], the transition form factors are de- 
rived from the constituent quark three-point function, 
the vertices of which are the weak interaction coupling, 
7^(1—75), and two phenomenological Bethe-Salpeter am- 
plitudes for the £?( s ) and / (980) mesons. While the 
B s can be parametrized with the leptonic decay con- 
stant (known from lattice-QCD simulations), the lat- 
ter is more problematic since the //» is poorly deter- 
mined. In an attempt to formulate a suitable scalar 
Jo (980) vertex function, we constrained its parameters 
by means of experimental quasi two-body branching frac- 
tions, _D( S ) — > /o(980)P, P = tt,K. The advantage is 

that the F+°~* fS (q 2 ) and F B °^ fc > (q 2 ) form factors, 

(m P2 )\rs^(l- l5 )b\B s ( Pl )) = 

F?^ fS (<l 2 )(Pi +P2)»+F B ^ f «(q 2 )( Pl - P2 )^ (51) 
can be calculated for any physical time-like momentum 
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TABLE V: As in Table III but for the helicity h = 0. 



/bo [MeV] 


LOVP 


230 

LOVPH LOVPH + C 


260 

LOVPH LOVPH + C 


290 

LOVPH LOVPH + C* 


U.C ( \ 

a 2 ' (nit) 

U.C ( \ 

a 3 (mb) 

U.C ( \ 

a 5 ' (m b ) 
a"' c (m b ) 
a 9 (mb) 


54.51 -i80.86 
5.86 + i 2.69 

-7.17-i3.14 
0.09 + i 0.03 

-9.31 -i 0.14 


160.2 - i 132.6 161.0 - i 129.4 
1.11 + i 5.01 1.08 + i 4.87 

— 1.12 — i 6.10 -1.08 -i 5.92 
0.02 + i 0.07 0.02 + i 0.07 

-9.07 -i 0.25 -9.07 -i 0.24 


165.7 - i 132.7 170.6 - i 129.5 
0.87 + i5.02 0.65 + i 4.87 

-0.81 -i 6.11 -0.53 -i5. 92 
0.02 + i 0.07 0.02 + i 0.07 

-9.06 -i 0.25 -9.05 -i 0.24 


171.8 -i 131.2 180.6 -i 129.2 
0.60 + i 4.95 0.20 + i 4.86 

-0.46 -i 6.02 0.04 - i 5.90 
0.02 + i 0.07 0.01 + i 0.06 

-9.05 -i 0.25 -9.03 -i 0.24 



transfer q 2 = (pi — P2) 2 ■ The superscript s is a re- 
minder that the transition is to the ss component of 
the scalar meson and p\ and pi are the B° and /o(980) 
four-momenta, respectively. We do stress that the B s — > 
/o(980) form factor calculated by us in Ref. [10] does not 
assume a pure ss state of the / (980). Instead, it was 
treated as a mixture of strange and non-strange qq com- 
ponents related by a mixing angle which also yields the 

related form factor F 01 (q ). This angle was deter- 
mined with experimental constraints [10] and the overall 
normalization of the transition form factor receives con- 
tributions from both states. 

The form factors F±(q 2 ) (we suppress the flavor super- 
scripts) are related to the set of vector and scalar form 
factors as, 



i^ 2 ) = F + {q 2 ) , 
F (q 2 ) = F+(q 2 ) + 



in: 



m 



(52) 

-F-(q 2 ) ■ (53) 



J! 

JO 



The form factor Fi(q 2 ) we obtain in both the dispersion 
relation and covariant light front dynamics approaches 
agree at the maximum recoil point q 2 = 0. At large 
four-momentum transfer, specifically for q 2 = Tny,^ ~ 

10 GeV 2 , our model predictions differ significantly which 
is also known to occur for B — > ir transition form fac- 
tors [32]. This is not surprising, as for large momentum 
transfers the final-state meson is less energetic and the 
soft physics of the bound states becomes more relevant. 
Since the models differ in their parametrization of the 
bound-state wave functions, it is clear that their inaccu- 
racies are revealed in the form-factor predictions at large 
q 2 . In Ref. [30], we deduce from the author's extrapola- 



0.3, which is 



tion parametrization that F 1 s (mj/^) 
compatible with our dispersion-relation prediction ~ 0.4 
within the errors. In Section VI, we will account for this 
rather large window of values and plot the ratio TZ^/^ as 

a function of F^'^^ {m 2 ,^). 



C. Model parameters 

The hard scattering contributions involve endpoint di- 
vergences, which we choose to parametrize by, 

X H = [1 + Ph cxp(i<t> H )) In -j^- . (54) 

In case of a possible annihilation or "other" contribution 
we simply write, 



X c = p c exp(i<f> c ) 



(55) 



which introduces four parameters, < pc.H an d < 



We assume that X, 



h=0 _ 



X 



<t>c,H < 360°. w 
X^r^ 1 = X C H , as the vector <fi an d scalar / (980) 
mesons have similar masses and we consider the ss com- 
ponent only. The hard scattering corrections are ex- 
pected to be of the order of m B o/\ h in Eq. (54), with 
A/j = 0.5 GeV. The parameters Pc,h and 4>C,h are cho- 
sen so as to reproduce the experimental data discussed 
in Section VI. We insert their values in the B° s — > J/V'/o 
decay amplitude (4) and then predict the branching ratio 
B(B° S -> foJ/,/,). 



VI. RESULTS AND EXPERIMENTAL DATA 

In the B® — > (f>J/tp decay, one can define five observ- 
ables: a longitudinal, parallel and perpendicular polar- 
ization fraction, and f±, respectively, 



fk 



|^ L | 2 + U|i| 2 + |^ ± 



- 2 , k = L,\\,± 



as well as two relative phases, cj)\\ and (f>j_, 



4>k = arg 



, k=\\,± , 



where we have abbreviated, Al = - / ^b°-7^/V' 



r A (h=+i) 



(56) 



(57) 

- ^11 = 
-i) 

4>J/4> 
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TABLE VI: Short-distance amplitudes, a q n (m b ) x 10 3 , for B° J/ipf (980) as a function of the f B o decay constant with 
f f s = 380 MeV and Ff°^ fS (m 2 J/4> ) = 0.4. See caption in Table III for the definition of LOVP, LOVPH and LOVPH + C 
amplitudes. 



/so [MeV] 


LOVP 


230 

LOVPH LOVPH + C 


260 

LOVPH LOVPH + C 


290 

LOVPH LOVPH + C 


U.C ( \ 

a 2 ' (m b ) 

U.C { \ 

a 3 ' (m b ) 
ar(m b ) 

U.C / \ 

a 7 ' (m b ) 

U.C / \ 

a g ' (mi,) 


11.61 - i80.86 
7.29 + i 2.69 
-7.17-i3.14 

0.09 + 0.03 
-9.38 -i 0.14 


-42.40 - i 255.5 -33.35 - i 224.3 
9.71 +i 10.53 9.30 + i 9.13 

-10.25 -i 13.12 -9.74 -i 11.34 
0.13 + i0.15 0.12 -M 0.13 
-9.51-i0.53 -9.49-i0.46 


-66.51 - i 234.1 -51.82 - i 224.4 
10.80 + i 9.57 10.13 + i 9.13 

-11.63 - i 11.90 -10.79 - i 11.35 
0.14 + i 0.14 0.14 + i0.13 
-9.56 -i 0.48 -9.53 -i 0.46 


-95.23 - i 229.5 -69.17 - i 223.8 
12.08 + i 9.36 10.91 + i 9.10 

-13.27 - i 11.64 -11.78 - i 11.31 
0.16 + i 0.13 0.15 + i0.13 
-9.63 -i 0.47 -9.57 -i 0.46 



TABLE VII: The phenomenological contributions ( h x 10 3 for h = 0,-1, +1, Eq. (24), to the B° -> J/ip<p amplitude as a 
function of the f B o decay constant with //« = 380 MeV. 



/bo [MeV] 


230 


260 


290 




-18.11 + i 1.98 


-28.04 + i 1.89 


-37.19 + i 2.63 




-129.26 + i 14.12 


-200.12 + i 13.46 


-265.41 + i 18.77 


c h=+1 


-15.25 + i 1.67 


-23.61 + i 1.59 


-31.31 + i 2.21 



The CP average is defined in terms of the polarization 
fractions, fk, 



A k — 



/ s £ s 

k ~ / fc 

J fc ^ J fc 



(58) 



Similarly, for B° -> / (980)J/V>, the CP average is de- 
fined as, 



iCP 



ff(p s ° -> /qJ/vq - g(P s ° -> /o J/VQ 

-+ /o J/V) + B(B0 -+ / J/V) 



(59) 



We use data from CDF and D0 for the B° -> <j>J/ip 
decay, whereas there is no available data on the channel 
B® — > foJ/ip. Our data compilation consists of the D0 
values for the amplitudes, |^ L | 2 = 0.555± 0.027± 0.006, 



0.244 ± 0.032 ± 0.014 and the relative phase 



0|| = 2.72± ;^ rad [7]. The CDF values [33] are compat- 
ible, |^ L | 2 = 0.530 ± 0.021 ± 0.007 and | 2 = 0.230 ± 
0.027±0.009, and the PDG data book quotes the branch- 
ing fraction, B(B° S -> J/ij}<j>) = (9.3 ± 3.3) x 10~ 4 [18]. 

The ratio 7£/ m has been argued [8] to be of the order 
0.2 — 0.3, based on the knowledge of the experimental 



ratio of decay rates [34], 

T(D+ -> / 7T + -> K+K-n- 
T(D+^ (t>TT+^ K+K-tt- 



0.28 ±0.12, 



(60) 



and an estimate of the scmileptonic, integrated branching 
fraction ratio B(D+ -> foe + v,f -> Tr + TT~)/B(D+ -> 
<£e + z/,<£ -> K+K~) = (13 ± 4)% from CLEO [35]. The 
ratio 7^/ o /0 was reassessed in terms of the differential 
decay ratio [36], 



dT 



n 



(D+ -> f e+vj a -> 7T+7T )| 
= 0.42 ±0.11. 



-> K+K- 



\q 2 =0 



(61) 



If we combine the above three experimental estimates, 
we propose a window of 0.2 < r Rf /^ > < 0.5 for the ratio 
based on D s decays. 

With the experimental data listed under Eq. (59) as 
constraint, we find optimal values for Xq and X^. In 
principle, we deal with a system of four coupled non- 
linear equations for \Al\ , \A\\ | , <fi\\ and B{B° S — > (j>J/tf>) 
and four variables, which does not put tight constraints 
on the phenomenological part of our B° s — > J/ip(j) am- 
plitude. When solving numerically we find, depend- 
ing on the fso values, two solutions among which only 
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TABLE VIII: Same as Table VII but for the B° -> J/^/o(980) amplitude. 



/so [MeV] 


230 


260 


290 


C 


-44.08 + i4.81 


-68.25 + i 4.59 


-90.51 + i 6.40 



TABLE IX: Prediction for the B° -> J/^/o observables for the different amplitudes LOVP, LOVPH and LOVPH + C along 
with experimental analysis data of the B° — > J/ip4> decay. Here central values, f B o = 260 MeV and //« = 380 MeV, and the 
transition form factor Ff 3 ^^° (q 2 = m 2 ^) = 0.4 are used. The values in the second column are predictions. Those of the third 
column include the hard scattering corrections with the endpoint parametrization pu = 1.85 ± 0.07 and <f>H = 255.9° ± 24.6°. 
The fourth column corresponds to the reproduction of the data with the parameters pu, 4>h, pc and <j>c displayed in the second 
line of Table II. 





LOVP 


LOVPH 


LOVPH + C 


Experimental 




(Prediction) 


(2 parameters) 


(4 parameters) 


data 




0.172 


0.554 


0.555 


0.555 ± 0.033 [7] 


\M 2 


0.404 


0.219 


0.244 


0.244 ± 0.046 [7] 


h (rad) 


-0.221 


2.13 


2.72 


2.72 ± 1.38 [7] 


B(B° -> J/ip<t>) 


0.00075 


0.00115 


0.00093 


0.00093 ± 0.00033 [18] 


B{B» -> J/^/ ) 


0.00020 


0.00047 


0.00050 




A C p{B» -> J/^/o) 


-0.00013 


-0.0013 


-0.0011 






0.28 


0.42 


0.55 





one yields a reasonable value for the branching fraction 
B(B S — > foJ/i/j) not too different from that in a naive 
quark model. We list the parameters pc,H and <j>c,H 
independent of ff* for three values of f B o in Table II, 
from which it is plain that the uncertainties on the mag- 
nitude of the modulus pc as well as the phase 4>c are 
substantial. The experimental errors on the observables 
are clearly not constraining enough. Yet, we observe that 
the variations of Xc and Xh are smooth as a function 
of the decay constant /bo . 

Likewise, we present numerical values for a^ h (mb) for 
the three helicities in B® — > J/ip(j) in Tables III, IV and 
V and for — > J/ipfo in Table VI as functions of /bo to 
illustrate one facet of the hadronic uncertainty. In these 
tables, we list the decomposition of a% h (rrib) for each 
value of /bo; in the first column, the values of a^ h (mb) 
are for the calculated leading order (LO), vertex (V) and 
penguin (P) amplitudes only. These are independent of 
/bo and correspond to the predictions in Figure 1. Next, 
the a^ h (mb) that contain the LO, V, P and the hard- 
scattering {H) amplitudes, where only pn and 4>h are 
fitted to reproduce the B® — » (f>J/ip observables while 
X c = 0. For f B o = 260 MeV one obtains p H = 1.85 ± 
0.07 and <j>H = 255.9 ± 24.6. These values are not very 



different from those given in the second line of Table II. 
This case corresponds to Figure 2. At last, denoted by 
LOVPH + C, we give the values for a% h (mb) for the case 
that the amplitudes are included, which corresponds 
to the pc,H and (j>c,H values in Table II and to Figure 3. 
We remind that the dependence on f B o enters the short- 
distance coefficients via the hard-scattering contribution 
H^{M\J/ip) in Eq. (22) and that the phenomenological 
amplitudes, Xh and Xc, are in competition with each 
other. Therefore, the hard scattering contributions to 
a% h (rrib) in LOVPH are slightly different than those to 
LOVPH + C. 

The largest values observed in the leading amplitude, 
a2 ,c {rrib), are for h = 0. We also remark there is 
no variation as a function of f Bs in Table IV since 
H% = ~ 1 (MiJ/'ip) = 0. Moreover, penguin contractions 
only contribute to a^' h ~ +1 (mb) and a|' ?l_+1 (?7ifc) in the 
Bg — > <j>J/ip amplitudes, while there are no penguin 
terms in B® — > f J/ip. Altogether, the penguin contri- 
butions are very small. We note that the contribution of 
the phenomenological amplitudes, (Tables VII and 
VIII), is small, about 6 — 7% of the h = 0, +1 amplitudes 
in B s — > (f>J/ip and 2% of the B s — > foJ/i/j amplitude, yet 
dominant in the h = — 1 amplitude devoid of penguin and 
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hard scattering corrections. Thus, any contribution from 
new physics, and to less an extent annihilation topolo- 
gies, should occur in the h = — 1 helicity amplitude. 

When including all the contributions (LOVPH + 
C), we qualitatively verify the hierarchy relation, 



]A (h=o) , 

l-^BO_j.0j/^| 

J/ip(j) and \A 



> 



(h=0) 



> \A 



B°- 



> \A 



(h=-i) 



-i) 
><t>j/i>\ 



in B { 1 



B°- 



, > U (/l=+1) I 
<j>J/il>\ ^ ^B°^<pJ/ip\ 



111 



the CP conjugate decay B° s — > J/ipcf). These hierarchy 
relations are also reproduced for the amplitudes when 
they include, besides tree contributions, vertex, penguin 
and hard-scattering corrections. 

Having determined numerical values for Xjj and Xq, 
we can calculate the B® — > fo J /ip amplitude and obtain 
the associated branching fraction and CP asymmetry. 
We do so for the central values of Jb'> — 260 MeV and 
ffs = 380 MeV discussed in Section V. For a transition 



form factor F x 
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) = 0.4 and for the differ- 



ent amplitudes LOVP, LOVPH and LOVPH + C de- 
fined above, those obscrvables are displayed in Table IX 
together with a comparison of the B° J/ip(j) results 
with the corresponding available experimental analysis 
values. Furthermore, we obtain for a transition form fac- 



tor F? 



(<1 2 



m 



) = 0.2: 



B{B S — > foJ/ip) = 3.80 x 10~ 4 , 
A C p(B s -> / J/V>) = -0.0005 , 
B{B S -> fij/ip) = 9.30 x 10~ 4 , 



n 



fo/ 



0.42 ; 



forFf^V = m^)=0.3, 



B(B S -> / J/V) - 4.37 x 10- 4 , 

A C p(B s -> / J/V) = -0.0008 , 

£(B S -> = 9.30 x 10~ 4 , 

^/o/0 = 0.48 ; 

and for if^V = = °- 5 > 

£(5 S / J/V) = 5.7 x 10- 4 , 

A C p(S s -+ / J/V) = -0.0013 , 

B{B S -> </.J/V) = 9.30 x 10~ 4 , 



K 



fo/ 



0.63 , 



and finally, the CP asymmetries in B s — > are, 



-3 



A£ P (B S -> </>J/V>) = -1.66 x 10 
Al P (B s -> </>J/V>) = 1.99 x 10" 
^cp( b s ~> = 2.15 x 10" 



Our prediction for the time-integrated asymmetry 
Acp(B s — > foJ/ii) is about one order of magnitude 
smaller than the Standard Model value, — 2/3 s = —0.036. 



We remark that the above numerical values for this CP 
asymmetry have to be interpreted with care — we choose 
the parameters of the full QCDF amplitude in Table II 
such that the experimental B° s — > J/ij)(j> observables are 
reproduced. In doing so, we may deliberately include 
"new physics" effects with just the Standard Model am- 
plitude, in particular via the additional amplitudes 
Moreover, we use the same end-point parameterization, 
Xh, in both decay channels since the B® — > J/ijifo 
branching ratio is not experimentally known. This ap- 
proach seems reasonable, as the physics buried in these 
infrared divergences must be similar in both decays. 
It could also lead to an overestimation of the hard- 
scattering contributions to B® — > J/V>/o as well as of 
A CP (B S -+ / J/V). 

We illustrate the variation of the ratio, 1Zf /^, by 
taking into account the uncertainties in the decay con- 
stants fgo and ffs as well as those in the decay 
rates, B(/ o (980) — > Tr+Tr") = 0.50tg^ [10, 36] and 
B{4> -> K+K~) = 0.489 ± 0.005 [18]. ' The results arc 
displayed in Figures 1-3. 

In Figure 1, 7£/ o /0 is plotted as a function of 

B°— >f 3 

F 1 3 (nij^) where only the tree amplitude along with 
vertex and penguin corrections are included in both am- 
plitudes, Ago^^j/^, and Asa^foJ/ip- The ratio is plot- 
ted with the corresponding envelope of 7£/ o /0 due to the 
uncertainty on the decay rates. In Figure 2, we aug- 
ment this amplitude by hard-scattering contributions, 
that is the full QCDF amplitude given in Eq. (22). Fi- 
nally, in Figure 3, 7£/ o /0 is plotted as a function of 



) including hard-scattering corrections and 



possible "other" contributions, Although the afore- 
mentioned uncertainties are depicted in all figures, we 
stress that those on the decay constants fgo and ff , 
where they apply, have more impact on the 72-/ o /0 band 
than the / (980) and <j> decay rate incertitudes. The 
spreading of the curves representing TZf /^, as a function 



of F- 



B^fS 



(m 2 . 



yi'jijj J j is larger with respect to the variation 
in fgo than in ffs . This points to the necessity of hav- 
ing an improved experimental determination of /go . The 
upper limit of the envelope is reached only for the largest 
values of fso and //» considered here. 

Figure 3 shows that our central-value predictions of 
72-/ o /0, in absence of any phenomenological contributions, 
are within the estimate by Stone and Zhang [8] for most 

values of the form factor F 1 3 {rrij,^). However, when 
the additional amplitudes, £, are accounted for in the 
decay amplitudes of Eqs. (2) and (4), the ratio 71^/^ 
exhibits three striking features: 

• Additional amplitudes, (, can play a major role 
due to their large contributions to both the numer- 
ator and denominator of the ratio 7£/ o /0, as seen 
from the comparison of Figures 2 and 3. 

• The predicted TLj Q j^, band overlaps well with the 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



FIG. 1: The ratio TZf /tj> as a function of the transition form factor F 1 (mju). Only tree, vertex and penguin contributions 
(LOVP), independent of the decay constants f B o and // , are included in the decay amplitudes. The dotted line corresponds 
to the central value of this ratio while the area between the two solid lines gives its envelope due to the uncertainties on the 
decay rates /o(980) — > tv + tt~ [10, 36] and <\> — > K + K~ [18]. The two horizontal dash-dotted lines delimit the (shaded) area 
between the experimental predictions found in Refs. [8] and [36]. 




0.1 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

TT^rt/o. 2 2 x 

r i (q =m m ) 

FIG. 2: The ratio 72. / /^ as a function of the transition form factor Ff s ^^ n {m 2 j^) where now the tree, vertex, penguin, and 
hard-scattering contributions (LOV PH) are included. The area between the two dashed lines gives the envelope of this ratio 
when taking into account uncertainties on the decay constants (f B o = 260 ± 30 MeV and // = 380 ± 40 MeV) while the 
solid lines include in addition the uncertainties on the decay rates /o(980) — ¥ n + n~ [10, 36] and <f> — > K + K~ [18]. The single 
dotted line is our prediction for the central values of the decay constants, f B o = 260 MeV and ff — 380 MeV. The horizontal 
dash-dotted lines correspond to the experimental predictions of Refs. [8] and [36]. 



estimates of Refs. [8] and [36] for Fj "^ to (m 2 J/ifj ) < 
0.4; beyond, our predictions are much larger, which 



may indicate a larger pollution due to /o(980) — > 
K + K~ if contributions from other than the Stan- 
dard Model were present. 
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• The uncertainties on the /o(980) and 4> decay 
rates increase the width of the band considerably, 
though the main uncertainty stems from the decay 
constants Jbo and // . 

Let us remind that the decay constant ffs only en- 
ters the hard-scattering and additional phenomenolog- 
ical contributions (C) to the decay amplitude B° — > 
/o(980) J/ip- If these are turned off, as in Figure 1, 
the ratio Rf /<p is still significantly above 10% for re- 
alistic transition-form factor values. That said, for 
practical purposes we decide to only consider the more 
recently obtained decay constants in Ref. [22] and 
choose three values within the given errors, //= = 
340, 380, 420 MeV. The value 180 MeV [25] yields too low 
branching fractions in other decays, for example B — > 
/o(980)7r,/o(980)p,/o(980)i<:W. Nevertheless, since we 
fix the hard-scattering parameters, pu and (f>H, only 
via the decay — > J/ipcf) and ff>> enters the numera- 
tor in 7?-/ o /0 linearly, the modification is straightforward: 
f fs = 180 MeV is about half the value //. = 380 MeV, 
therefore the central value of 7\L/ O in Figure 2 decreases 

from 0.42 to 0.19 (for F^ f " = 0.4). This is still 
within the limits predicted by the experimental esti- 
mates, 0.2 < Hf /j> < 0.5, and implies an 5-wave pol- 
lution. 

We infer from our numerical results that S'-wave kaons 
or pions under the <j> peak in B® — > J/ip<fi are very likely to 
originate from the similar decay — > J/if)fo- Therefore, 
the extraction of the mixing phase, — 2/3 s , from B® — > 
J/ip<t> may well be biased by this S- wave effect which 
should be taken into account in experimental analyses. 
In our interpretation of the full QCDF amplitude, we 
not only confirm the influence of S'-wave contamination 



as advocated in Refs. [8] and [36] but also find that its 
effect could be sizable. 



VII. CONCLUSIVE OUTLOOK 

The "phase" of B° - B° s mixing, -2/3 s , is thought to 
be best measured in the golden decay, B® J/ip'P: an d 
provides an opportune place to investigate physics be- 
yond the Standard Model. Several models have been pro- 
posed to explain the apparent discrepancy of the Stan- 
dard Model prediction for — 2j3 s with recent experiments, 
in particular exploring the impact of heavy, as of yet 
undiscovered particles on CP violation in weak B- meson 
decays. A general analysis of possible new physics ef- 
fects in the case of B° — B® mixing was recently given 
by Chiang et al. [21]. In there, the authors investi- 
gate several beyond Standard Model variations of the 
B s — > J/V"^ decay, such as ^-mediated Flavor Chang- 
ing Neutral Currents (FCNC), two Higgs doublets and 
SUSY, and find that new physics contributions may only 
modestly contribute to the mixing phase. However, it 
is also concluded, somewhat prematurely, that the CDF 
and D0 results are clear signs of new physics. 

In the present paper, we have taken a different path 
and studied the contamination of final state S-waves 
kaons in the B® —> J/ipcf) channel by those originating 
from the / (980) in the very similar B° s -> J/i/>/ (980) 
decay. We find that this effect is strong enough already 
for amplitudes including leading order, vertex and pen- 
guin corrections to create a real bias in the determination 
of -2)8.. 

Of course, we are aware that the phenomenological 
endpoint parametrization of a s corrections in the am- 
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plitudes H n (M\J /iji) and H%(MiJ/i/j) can cloud possi- 
ble new physics contributions alongside the C^ 1 ) contri- 
butions. In this case, we suppose that any new effects 
should be of comparable magnitude in — > J/ijxj) and 
B° s — > J/-0/ o (98O). Therefore, the S- wave contamination 
would be on the upper side of the estimate we propound 
and future analyses of the mixing angle in B s decays 
should be concerned with this effect. 
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